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ABSTRACT 

Alfven waves may be generated via mode conversion from fast magneto-acoustic waves near their 
reflection level in the solar atmosphere, with implications both for coronal oscillations and for active 
region helioseismology. In active regions this reflection typically occurs high enough that the Alfven 
speed a greatly exceeds the sound speed c, well above the a — c level where the fast and slow modes 
interact. In order to focus on the fundamental characteristics of fast/Alfven conversion, stripped of 
unnecessary detail, it is therefore useful to freeze out the slow mode by adopting the gravitationally 
stratified cold MHD model c ^ 0. This provides a benchmark for fast-to- Alfven mode conversion in 
more complex atmospheres. Assuming a uniform inclined magnetic flcld and an exponential Alfven 
speed profile with density scale height h, the Alfven conversion coefficient depends on three vari- 
ables only: the dimcnsionless transverse-to-the-stratification wavenumber k. = kh, the magnetic field 
inclination from the stratification direction 9, and the polarization angle (j) of the wavevector rela- 
tive to the plane containing the stratification and magnetic field directions. We present an extensive 
exploration of mode conversion in this parameter space and conclude that near-total conversion to 
outward-propagating Alfven waves typically occurs for small and large (f> (80°-90°), though it is 
absent entirely when 9 is exactly zero (vertical field). For wavcnumbcrs of helioseismic interest, the 
conversion region is broad enough to encompass the whole chromosphere. 
Subject headings: Sun, oscillations; sunspots; Magnetohydrodynamics (MHD) 
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1. INTRODUCTION 

Magnetohydrodynamic (MHD) linear mode conversion is an important process in solar active regions and in the 

overlying atmosphere Fast-to-slow conversion is implicated in the absorption of p-modes by sunspots (Cally et al. 
1994; Crouch ct al. 2005) and is well-understood in terms of local analysis around the Alfven/acoustic equipartition 
level a = c where the soimd speed c and Alfven speed a coincide (Schunker & Cally 2006). 

To be specific, our primary target in this paper relates to mode coupling in sunspots between the well-known p- 
mode seismic wave field of the solar interior and oscillations of various types in the overlying atmosphere. Sunspot 
seismology is one of the most difficult issues confronting helioseismology today (Moradi et al. 2010). Early attempts 
at seismic inversion using time-distance helioseismology interpreted phase shifts in terms of temperature variations in 
the first few Mm below the surface, and obtained a two-layer thermal structure that is at odds with other inversions 
(see Figure 19 of Moradi et al. 2010). Clearly, the near surface layers of sunspots are dominated by magnetic fields, 
with the plasma (3 typically passing through unity a few hundred km below the surface in umbrae, and at about the 
surface in penumbrae. This has the effect of both mandating fast-to-slow mode conversion where the sound and Alfven 
speeds are equal (Schunker & Cally 2006; Cally 2007; Cameron, Gizon, & Duvall 2008), and radically changing the 
phase of the fast wave that emerges through the surface and then reflects back downward to rejoin the subsurface 
helioseismic field (see Cally 2007, Figure 3, and Cally 2009). The combination of these two effects must confound any 
seismic inversion effort that attributes phase anomalies to sound speed perturbations alone. Conversion between the 
fast and slow magnetoacoustic waves may occur in two dimensions (2D), where magnetic field and the wavevector lie 
in a vertical plane {x-z say), and also in 3D. 

A further insufficiently modelled effect associated with this picture is that of fast-to-Alfven conversion, occurring 
near the fast wave reflection point when the waves are directed at an angle to the magnetic plane (3D). Typically, 
this occurs a few hundred km above the a = c level. Conversion was quantified for uniform inclined field in a specific 
simple solar atmospheric model by Cally & Goossens (2008) and confirmed in simulations by Khomenko & Cally 
(2011). The process has the potential to both remove energy from the reflecting fast wave and alter its phase before 
it re-enters the interior, with obvious implications for helioseismic inference, and for our interpretation of observations 
of waves in sunspot atmospheres. This explains our concentration on fast-to-Alfven conversion. However, the reverse 
Alfven-to-fast process is governed by the same considerations. It is well-known that the transition region can act as a 
powerful reflector of Alfven waves (Hollweg 1981; Cranmer & van Ballegooijen 2005), so escaping Alfven waves may 
in fact be partially 'reabsorbed' by the fast wave field after such reflection, though presumably with radically different 
phase. This would further complicate the seismology of sunspots. Alfven reffection and reabsorption is beyond the 
scope of the present study though. 

In light of this context, it is appropriate to concentrate on the vertical rather than horizontal variations in Alfven 
speed in the region of interest. The Alfven scale height in a low sunspot atmosphere may be of the order of 100- 
200 km typically, whereas horizontal variations may take place over distances of order 10 Mm. In the interests 
of mathematical tractability, we therefore ignore the latter. However, full 3D numerical simulations in a realistic 
spreading sunspot magnetic fleld have very clearly conflrmed the broad conclusions to be presented here (Khomenko 
& Cally, in preparation). 

Irrespective of these considerations, the basic process of fast-to-Alfven conversion is of fundamental interest in MHD 

wave theory, and is pursued here in that spirit. 

The full details of speciflc atmospheric models may be only weakly relevant to the strength of fast-to- Alfven mode 
conversion, possibly implicated in the apparent ubiquity of Alfven waves^ in the solar corona, as observed using 
the Solar Optical Telescope (SOT) aboard Hinode (De Pontieu et al. 2007) and with the Coronal Multi-Channel 
Polarimeter (CoMP) at the National Solar Observatory, New Mexico (Tomczyk et al. 2007). Figure 1 displays a 
typical example where fast-to-slow conversion occurs near the a = c equipartition depth but the fast and Alfven modes 
come to near- coincidence considerably higher, where c <C a. This suggests that we seek to simplify the analysis of 
fast-to- Alfven mode conversion by 'freezing out' the slow wave. This can be done by examining a cold MHD plasma 
c = in which the uncoupled fast wave dispersion relation becomes simply = a^|kp whilst the Alfven wave remains 
uj^ = a^fcj^, where oj and k are the frequency and wavevector as usual, and the subscript || indicates a component in the 
direction of the magnetic fleld . The cold MHD approximation corresponds to the ,6 = limit, where /? is the ratio of 
gas to magnetic pressure. We briefly discuss the applicability of the cold-plasma approximation to solar atmospheric 
mode conversion in Section 6 by comparing with some relevant warm plasma results. 

One advantage of the /3 — approximation is that we may ignore explicit gravity and just retain the density 
stratiflcation it produces. In large scale magnetic structures such as sunspots the density scale height is typically much 
smaller than the scale of variation of magnetic fleld Bq. We may assume Bq is uniform and therefore that the variation 
in the Alfven speed a = Bq/^JJTp is entirely due to the density stratification. Our aim in this paper is to quantify 
fast-to- Alfven mode conversion for this simplifled cold plasma model in the expectation that it will apply broadly to 
more complex structures and atmospheres. 

This model has already been partially explored by Cally & Andries (2010) and Hanson & Cally (2011), who examined 
the situation where the magnetic fleld is oriented perpendicular to the direction of stratiflcation x and the density 
decreases exponentially with scale distance h, i.e., oc exp[a;//i]. In the sunspot model this would be horizontal 

^ There has been some controversy in the literature about whether these are truly Alfven waves, or instead kink waves (Van Doorsselaere 
et al. 2008). Kink waves occur readily in media with transverse variations in the Alfven speed, for example in atmospheres consisting of 
isolated or packed flux tubes, and propagate at a phase speed which is a weighted average of the internal and external Alfven speeds. The 
models we address here have no such transverse structure, and so are undoubtedly Alfven waves. They may well be expected to couple to 
kink waves higher in the atmosphere though where coronal loop structure dominates. 
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Figure 1. Typical z-k^ propagation diagram (wtiere 2 is height) for a 5 mHz wave in the model of Cally and Goossens (2008) where the 
magnetic field is inclined at angle 8 = 25° from the vertical and the wavevector is oriented </> = 30° out of the vertical magnetic plane. 
The inner lobe represents the fast wave and the outer branches correspond to the slow wave. The intermediate curves are the Alfven wave. 
The vertical line represents the equipartition depth at which a = c. Note that the fast wave reflects at about 200 km whilst the slow and 
Alfven waves extend indefinitely upwards. The slow wave does reflect at lower frequencies below the ramp-reduced acoustic cutoff oJc cos 9, 
though that is not relevant to our discussion here. 



magnetic field and may be relevant to penumbra and canopy. It is a singular case in which resonant absorption 
occurs at the Alfvenic critical level where = a^/c^. Nevertheless, Cally & Andries succeeded in interpreting this 
absorption as a mode conversion, thereby clarifying the relationship between these two related processes. With y and 
z wavenumbers ky and fixed, they tabulated and plotted the fast-to- Alfven conversion coefficient ^(cr, (/>), where 
tan0 = ky/kz denotes the wave polarization, a = k^/^ sin^ </>, h is the uniform density scale length, i.e., oc exp[a;//i], 
and K, = {ky + klY^'^h = (k^ + k^)^/^ is the dimensionless transverse wavenumber. Remarkably, ^ is limited by a 
maximum value = 0.4937 attained at ct = cjo = 0.4644 in the limit (p ^ {] (see Figure 2a). Note though that R = 
in the 2D case </> = 0; Alfven conversion requires a 3D component Ky in the wavevector. The seemingly contradictory 
result that absorption is maximal in that limit in fact applies only if a is kept fixed at ctq and hence that k ^ 00 
simultaneously. 

Our task here is to generalize this result to arbitrary magnetic field orientation 6 from the x-direction. Without loss 

of generality we may assume that Bo lies in the x-z plane but that the incident wave may cut across it (ky ^ 0). This 
is more straightforward than the perpendicular field case since there is no resonant singularity in the equations. 

2. EQUATIONS 

Consider a cold MHD plasma with monotonic increasing Alfven speed a{x) (due to monotonic decreasing density) and 
uniform magnetic field Bo = -Bo (cos 0, sin 6') in cartesian coordinates {x,y,z), with 0° ^ 9 < 90°. The orientation 
of the stratification is arbitrarily chosen to be in the x rather than the z direction for consistency with Cally & 
Andries (2010). In the linear approximation the plasma displacement ^(a;, y, z, t) = ^{x) exp[i{kyy + k^z — wt)] where 
^{x) = ^x^x + ^yGy + ^zGz = Cj-Gj. + ^yGy then satisfies 

(5||+^)^ = -^pX, (1) 

generalizing equation (20) of Cally & Andries (2010). Here x = is the dilatation, the subscript denotes the 
parallel direction (cos 6, 0, sin 9), 'j.' indicates the direction (sin 9, 0, — cos 9) perpendicular to the field in the x-z plane, 
and 'p' refers to the component in the plane perpendicular to Bq. Since there is no restoring force in the parallel 
direction it follows that =0, and hence 4 = €p- 
The corresponding wave-energy (Poynting) flux is 

F = ^Re[EtxBi] = Folm [xt + it' 5||4)e,i] (2) 



4 



P.S. Cally and S.C. Hansen 



0,5 
0,4 
0,3 
0,2 
0.1 
0.0 



■ / 


^ . . . . 

\ 
— 










; / 

-nf 


\ 

5 










: 1 




i 

\ 

—\ 








1 

■1 




\ 








f 













0.0 



0.5 



1.0 



1.5 
cr 



2,0 



2.5 



3,0 




Figure 2. Left panel: Absorption coefficient J? as a function of cr = {kh)^f^ siv? (f> for oc e~^^^ and = 5°, 10°, . . . , 85° (thin curves, 
top to bottom) for the resonant case d = 90° (from Cally & Andries 2010). The heavy dashed curve represents the limit tf> 0° . The 
vertical line is at c = co = 0.4644. Right panel: Contour plot of A against </> and k = kh. The maximum in absorption occurs as k oo, 
— >■ along the curve k^/^ sin^ (j> = (delineated in white). 



where Fq = luBq/^, and Ei = — vxBq and Bi are the perturbed electric and magnetic fields respectively. It may be 
verified directly (with the aid of the product rule) that V-F = by contracting Equation (1) with ^* and taking the 
imaginary part. Since F is clearly independent of y and z it follows that the x-component is constant, as one would 
expect. This is distinct from the resonant case 9 — 90° where F^ is only piecewise constant, with a discontinuous drop 
to zero at the Alfven resonance. 



Eliminating x from (1) results in 



or in terms of x-derivatives only (for computational purposes). 



[cos Odx + ikz sin 0)"^ + ^ ~ kl ) 



-i kyd±^± , 

-iky (sin 6 dx — ikz cos 0)S^y 
-iky (sin 6 dx — ikz cos 0)^j_ 



(3) 



(4) 



In general, the fast and Alfven waves are intricately intertwined in these equations. Their separate identities are 
more clearly seen in Equation (1), where the left hand side exhibits the pure Alfven operator and the right hand side 
represents the fast wave (characterized by the dilatation x) &s a source term. The perturbation analysis of Section 5 
further expands on the coupling between the two wave types. 

We now specialize to the Alfven profile defined by uP^'b? ja^ — e^^/'', where h is the density scale length, and also 
define the dimensionless variables s = e~^^^, X = ^/h, Ky — kyh — Ksinc/), and Kz = k^h = Kcoscj). Note that fast 
mode reflection occurs (classically) where kx — 0, i.e., at — a^iky + kf). In dimensionless form this is s = k^. The 
Alfven wave is well-described by the eikonal approximation where s = k| 3> 1. 

Defining U ~ {X^,Xy, sAT^, sA^), where the prime denotes the s derivative, Equations (4) take the form 

sU' = AU (5) 

with 



I 

P Q 



where and I are the 2x2 zero and identity matrices respectively, 

^2 



-K COS sec sine/) k (tan^ + sin' 



K cos tl cos (p sm 
la I „;„2 



s sec 9 



(6) 



(7) 
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Figure 3. Same as Figure 2a, except for 6 = 85° (full curves). For comparison, the Alfven absorption curves of Figure 2 for the resonant 
case 6 = 90° are overlaid as dotted curves. The full curves are truncated at k = 8 for numerical reasons, which is progressively more 
restrictive as <p increases. 



and 

n — ■ f ^ sin6'sin^ \ , . 

^ ~^'^\secesm(j)taxie 2cos(j)tanej ' 

Note that we may write A = Aq -|- sAi, where Aq and Ai are constant matrices. 

Equation (5) may be solved numerically subject to the boundary conditions (i) the fast wave decays as x ^ +00 
(s — 0+); (ii) there is no incoming Alfven wave at x = +00; (iii) there is no incoming Alfven wave at a; = —00; and (iv) 
the incoming fast wave at a; = —00 carries unit .T-flux. To apply these conditions in practice we develop a Probenius 
expansion about s = and a WKB solution valid for large s in the Appendix. 

3. NUMERICAL RESULTS 

The resonant perpendicular field case = 90° discussed in Cally & Andrics (2010) is commonly thought to differ 
fundamentally from the non-resonant cases 0° < < 90° addressed here in that it is singular at the Alfven resonance 
w = akii , and the 'Alfven absorption' is generally perceived as a local resonant absorption rather than a simple mode 
conversion to a transmitted wave. (We have defined as the fraction of the incident fast wave energy flux that is 
ultimately converted to an Alfven wave. Fraction 1— R belongs to the reflected fast wave.) Nevertheless, Figure 3 for 
9 = 85° suggests that the two situations are not that different after all: the Alfven absorption/conversion coefficients 
arc near-identical. Clearly, despite the singular nature of the 9 = 90° 'resonant absorption' case, it is a continuous 
extension of the 9 < 90° 'mode conversion' cases. However, this is not to say that resonant absorption is illusory. 
See the discussion in Section 6 of Cally & Andries (2010) for further detail, in particular regarding the significance of 
trivial Fourier-transformable directions. 

Figures 4 and 5 show graphically how the forward Alfven conversion coefiicient (conversion to the rightward 
propagating Alfven wave) and reverse coefficient (conversion to the leftward Alfven wave) vary with transverse 
wavenumber k and wave polarization (f> for magnetic field inclination 9 = 10°, 20°, . . . , 80°. 

Clearly, is favoured at (j) < 90° and generally dominates on (/> > 90°. This is as expected, since maximal 
mode conversion is associated with an alignment of the phase velocities of the donor and recipient wave. If this occurs 
on the 'upstroke' of the fast wave's path then is stronger, but if it is on the 'downstroke', after reflection, then 
is the more significant. Since the approximate alignment occurs before the fast wave refiects for 0° < ^ < 90° and 
—90° < 4> < 90°, and after for cj) > 90°, the numerical results are plausible. The increasing symmetry between and 
as they weaken with increasing 9 is also consistent with this interpretation. 

Both vanish in various limits, specifically (i) k = 0; (ii) 9 = 0°; and (iii) </> = 0° and 180°. The contour figures 4 
and 5 indicate though that the drop to zero in the k ^ limit is very sharp for small 6. This is also seen in Table 1 
for at 61 = 10°. 

Figure 6 (top) shows a single frame from an animation for k = 1, 6 = 30°, (j) = 40° available as a supplement to this 
paper. It depicts the fast wave field through the green-yellow background shading representing x- Several magnetic 
field lines are overplotted, and wave back and forth in the animation. In the still frame shown here it is clear that to 
the right of the reflection point (a; = in this case) the field lines are nearly straight, because of the very long Alfvcnic 
wavelength, but tilted and moving with respect to each other. In contrast, the second supplemental animation (Figure 
6 bottom) for the same case but with cj) = 140° (exactly equivalent to reversing 9 but keeping (p unchanged) exhibits 
very little Alfvenic action, as we may surmise from Figure 4. 

Recognizing that wave direction </> is likely to be uniformly distributed in most circumstances. Figure 7 plots the 
0-averaged Alfven conversion coefiicient {Ji'^) against k for a range of field inclinations 9. 



4. FAILURE OF LOCAL ANALYSIS 



6 



P.S. Cally and S.C. Hansen 




Figure 4. Left column; The forward Alfvcn conversion coefficient j?+ as a function of 4> and k for 8 = 10°, 20°, 30°, and 40° (top to 
bottom) as labelled. Right column: The reverse Alfven conversion coefficient Sl~ for the same cases. In all cases, the contours are 0.1, 0.2, 
. . . , 0.9. The respective maxima in Ji+ are: 0.951 for = 10° (k = 0.056, </> = 88.0°); 0.814 for 9 = 20° (k = 0.107, if> = 82.9°); 0.647 for 
e = 30° (ft = 0.170, (f> = 73.9°); 0.505 for 9 = 40° (k = 1.29, </> = 41.0°). 



Mode conversion between fast and slow magnetoacoustic waves is amenable to local analysis (Schunker & Cally 2006) 
using the method of Tracy et al. (2003) or related WKB-based techniques since the conversion region is characterized 
by an avoided crossing of generic saddle point topology in the appropriate phase space. Tests against exact solutions 
in a gravitationally stratified isothermal atmosphere suggest very satisfactory accuracy when the gap between modes 
is small to moderate (Hansen & Cally 2009). As already indicated in Figure 1 though, fast /Alfven interactions 
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Figure 5. Same as Figure 4, but for 9 = 50°, 60°, 70°, and 80° (top to bottom). The maxima in J?"*" occur above ft = 5 in all cases. 



typically involve a long and definitely non-local conversion region which may extend for several scale heights and 
which does not display the required saddle structure. Figure 8 illustrates the phase structure for the cold plasma 
case at hand, again displaying the long distributed interaction between the fast and Alfven waves. The dispersion 
relation D = {s — \k\'^){s — Ky) = clearly exhibits the two apparently disjoint modes, which in s-k^ space are simply 

non-intersecting parabolae. Attempts at a local analysis around the point of closest approach of the two branches have 
produced unsatisfactory results. It does not appear that any such local analysis can adequately capture fast-to- Alfven 
mode conversion in this scenario. Indeed, surprisingly, we shall see in Section 5 that the bulk of the mode conversion 
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Figure 6. Top: A greyscale version of a frame from a supplemental animation for the case k = 1, 9 = 30°, </> = 40°. The shading 
represents x = V- ^, and hence the fast wave, which classically reflects at x = (vertical black line). A selection of magnetic field lines are 
overplotted. In equilibrium they are all inclined at 30° but when perturbed they oscillate up and down due to the Alfven wave. Beyond 
about X = though they remain almost straight due to the increasingly long Alfven wavelength. In the animations, the field line colours 
represent back-and-forth motion in the y direction. Bottom: Same, but for </< = 140°, exhibiting much less Alfven conversion. 



Table 1 

Forward Alfven conversion coefficient for a range of k and 4> at magnetic field inclination 6 = 10°. Similar tables for 
9 = 20° , 30° , . . . , 80° are provided as supplementary material. 
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Figure 7. averaged over all The various curves are for 6 = 10°, 20°, . . . , 80° (top to bottom). 




Figure 8. The x-k^ phase plane for k = 1, = 30°, <j> = 40°. The inner lobe represents the fast wave and the two outer wings are 
the Alfven wave. The grey shading corresponds to \kx\ < 1. Formally the eikonal approximation is valid only on \kx\ ^ 1, though the 
dispersion curves shown here are still instructive at small Kx , correctly indicating reflection of the fast wave at around a; = and continued 
propagation of the Alfven wave as x — > +oo. 



actually occurs around and beyond the fast wave reflection point rather than where the fast and Alfven loci approach 
most closely. 

5. PERTURBATION ANALYSIS 

Although the numerical analysis is essentially complete, it is of interest to develop a perturbation solution that 
illustrates the nature and locality of the fast-to- Alfven interaction at the core of the conversion process. We perturb 
about the purely fast case Ky = 0, — 0. Equation (1) is rewritten 



kI + e 



Ky^y . 



')?_L = -i Ky dj_^y 

(9^ + e-=")Cy = -i ny d^U 
The fully reflective fast wave with Kj, = and = n has solution 



_L0 



(9) 



(10) 



as a Bessel function of the first kind (standing wave), or alternatively in terms of leftward and rightward propagating 
Hankel functions. The incident 3;-component of wave energy flux associated with ^J|q is simply Fin = Fq/Att. 
The transverse displacement appears at first order in Ky. The correction ^j_2 to is of second order and will not 
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Figure 9. Perturbation theory Alfven conversion coefficient (full curve) given by Equation (15) as a quadratic function of Ky for 
K = 0.2, 6 = 30°. The exact numerical result is shown dashed. For this case B{oo) = —0.6539 — 1.257 «. 



be required. At first order 

iyi ~ -^e-*'=-^*^'^^sec2 0{>l(x)if^''(2e-"/2sec0) + S(a;)/f^''(2e-"/2sec0)} , (11) 

where 

/'QO 

A{x)= / e'^^^'^''^H^''\2e-^/^sece)d^^^o{X)dX , 

J X 

J — oo 

The limits on the integrals have been chosen so that there are no incoming Alfven waves at either end. For large x 
well beyond the interaction region, the upper limit of the B integral may be replaced by oo, whence 

Ci ^ -^sec2^B(oo)e-''^»^*^"^i?^'^(2e-*/2sec^) as a; ^ oo, (13) 

carrying x-component of flux 

2 

F+ ~ FoIm[Cyd\\^yCose] = Fq'^ |B(oo)|2 sec^ ^ . (14) 

Hence 

K^TT^sec^ 6*15(00)12 + 0(4)- (15) 

Figure 9 illustrates how well this quadratic result compares with the full numerical solution for small Ky. 

The interaction integral B{x) defined in Equation (12) provides a convenient picture of where fast-to- Alfven mode 
conversion occurs. The jump in \B\'^ determines the amoimt of mode conversion, and the argimient of B{oo) represents 
a phase shift. ^ Figure 10, for the case k = 0.2, 9 = 30°, shows that it is in fact far more spread out in x than might be 
expected from a dispersion diagram, with the major contribution occurring around and beyond the fast wave reflection 
point, where the eikonal approximation breaks down. This is consistent with the findings of Cally & Andries (2010) 
for the = 90° case. At higher k though, as one might expect, the growth in jSp is sharper and therefore the process 
more localized, but these higher wavemimbers are of lesser helioseismic interest. 

Figure 10b indicates that mode conversion is progressively more localized as k increases, as we might expect, though 
in all cases it predominantly occurs beyond the fast wave reflection point. However, since k < 0.2 typically for 
oscillations relevant to local helioseismology,^ it is clear that fast-to- Alfven conversion is spread out over many scale 
heights above the reflection point for these waves, easily encompassing the whole chromosphere. Although the present 
model does not contain a transition region (TR) or corona, we may surmise that the jump in temperature of order 100 
across the TR will sharply turn off this interaction. This is because the scale height h and hence k = kh are similarly 
increased by two orders of magnitude. 

^ The integrand of B is extremely oscillatory and hence there is a severe loss of precision in integrating it numerically. In calculating 
B(oo) it is advisable to deform the integration contour (—00,00) some distance ( < 27r to remain on the correct Riemann sheet) below 
the real line in the complex plane, where the oscillations are suppressed. This is done in calculating the B(oo) used in Figure 9 but, for 
purposes of illustration, Figure 10 uses integration on the real line. 

Assuming a density scale height of 150 km, the helioseimic degree i = 4640 k. 




(12) 
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Figure 10. Left: Interaction integral Ti{x) for k = 0.2, B = 30°. The real and imaginary parts are shown as full and dashed curves 
respectively. Right: k'^\B(x)\^ for k = 0.2 (full curve), k = 1 (dashed), and k = 5 (dotted). The vertical line indicates the position of fast 
wave reflection in the ft = 0.2 case. Reflection is at x = for k = 1 and x = —3.22 for k = 5. 




Figure 11. The lower panel of Figure 2 of Cally & Goossens (2008), showing the total magnetic (Alfven) flux (absolute, not relative to 
an incident flux) at the top of a simplified (warm) solar-like atmospheric model threaded by a 2 kG uniform magnetic field and topped by 
an isothermal slab above z = 0.5 Mm. The nine contours are equally spaced in flux. The bottom of the acoustic cavity hosting the waves 
in question is at 21 = —5 Mm, resulting in a horizontal wavenumber k = 1.37 Mm~^ for 5 mHz oscillations. This corresponds roughly to 
K = 0.2 in our dimensionless units. 

6. COMPARISON WITH /3 > MODEL 

The aim in this paper is to fully characterize fast-to- Alfven mode conversion in a cold MHD plasma with uniform 
inclined magnetic field and an exponential Alfven speed profile, neglecting slow-wave (acoustic) effects. However, it is 
instructive to compare our results with those of Cally & Goossens (2008) for a similarly configured warm plasma solar 
atmosphere to judge the extent to which the cold plasma model pertains to more realistic atmospheres. Figure 11 
shows the magnetic wave-energy flux emerging at the top of this model for a case roughly comparable to our k — 0.2. 
This is absolute flux (not relative as we have focussed on here) resulting from a driving plane at Zf, = —4 Mm and 
normalized by the total acoustic energy in Zf, < z < 0. Nevertheless, we might expect this figure to broadly coincide 
with the cold-plasma results. 

It is therefore surprising to view Figure 12, where strong Alfven conversion persists almost all the way down to 
= 0°, in stark contrast to the warm-plasma case where it is strongly suppressed for < 15°. Despite this, the 
behaviours at larger 9 are comfortingly similar, with strongest response around = 50°-70° in both cases. 

The reason for the disparity is made clear by the 5 mHz warm-plasma dispersion curves^ presented in Figure 13. 
At 6* = 40°, the fast- Alfven interaction is much as in the cold plasma model. The slow wave locus is sufficiently far 
away as to have no effect. However, aX 9 = 10° the slow locus impinges on the Alfven branch and actually causes it to 
turn over, suggesting refiection. The low-0 truncation in Figure 11 therefore makes sense; there will have been strong 
fast-to-Alfven conversion as in the cold plasma case, but it quickly reflects with probably some weak Alfven-to-slow 
(i.e., acoustic) conversion accompanying it. At lower frequencies, below ujcC0s9, the slow wave locus also turns over 
due to the acoustic cutoff effect. However, at higher frequencies (right panel) the generic cold plasma scenario is again 
apparent, even at small 9. 

7. CONCLUSION 

^ The dispersion function used here is that derived by Newington & Cally (2010), viz., 

V = uj'^uila^kl sin^ 9 sin^ + (aj2 _ o?kl) [w" - {o? -|- c^)ui'^k'^ + a'^c'^k^kf^ + c^N^k^ - {lj^ - a^k^ cos^ e)uil] , 

where ujc is the acoustic cutoff frequency, N is the Brunt- Vaisala frequency, a and c are the Alfven and sound speeds, fc = |k| is the total 
wavenumber, and fc^ is the horizontal wavenumber. The first (additive) term on the right hand side breaks the separability of the Alfven 
{i.e., u)^ — a^fc?) and magnetoacoustic (expression in square brackets) modes, but it vanishes in the cold plasma regime where ojc = 0. 
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Figure 12. .J?"*" against and 9 for fixed k = 0.2, roughly corresponding to Figure 11. As usual, tlie contours are 0.1, 0.2, 
that J?+ = on e = 0. 
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Figure 13. Warm-plasma dispersion curves corresponding to 6 = 40° (left) and 8 = 10° (middle) with (j) = 60° for the 5 mHz case of 
Figure 11. The inner lobes represent the fast wave, the outer wings the slow wave, and the intermediate branches the Alfven wave. The 
right panel corresponds to higher frequency (8 mHz), again with 6 = 10° and <j> = 60°. 



This paper addresses the fundamental issue of mode conversion between fast magnetoacoustic waves and Alfven 
waves in a cold stratified atmosphere. This is not to claim that the cold plasma approximation necessarily provides 
a full description of the solar atmosphere in conversion regions. It does though focus our attention solely on the 
conversion process. As indicated in Figure 13 warm plasma effects can further intrude to modify the ultimate Alfvenic 
transmission, and these must be taken into account in a full description of wave propagation through these regions. 
Nevertheless, they are distinct processes, and warrant individual attention. For example, at small field inclination there 
may be very significant Alfven conversion, but the Alfven waves may subsequently be reflected back downward. This 
complex array of distinct but spatially adjacent processes has implications for the wave field in the solar atmosphere 
overlying active regions, and also for the helioseismic wave field beneath them. 

Although we have exclusively focussed on fast-to- Alfven conversion, the reverse process may also be relevant in 
certain localized instances in the solar atmosphere. As a reversible process, conversion coefficients either way must be 
the same. 

For the cold plasma at least, reference to Figure 7 indicates that a randomly oriented incident fast wave field with 
K < 1 encountering a moderately inclined magnetic region, such as is found in sunspot umbrae, might be expected to 
lose of the order of 20% or more of its energy to Alfven waves. Given the strong </) dependence though, we may also 
conclude that such inclined field can act as a powerful directional filter on the waves entering the solar atmosphere 
through active regions, suggesting possible observational tests. 

In summary, the major insights we have gained from this study include: 

1. Conversion to Alfven waves is very sensitive to the direction (f> of the incident fast wave, as illustrated in Figures 
4 and 5. 

2. Unlike fast-to-slow conversion, the interaction between fast and Alfven waves is significantly spread across many 
scale heights. The major contribution to the conversion occurs around and beyond the fast wave reflection point 
rather than in the neighbourhood of the closest approach of the respective loci in phase space. 

3. For transverse wavenunibers k — kh < 0.2 typical of local helioseismic waves, the conversion region is sufficiently 
spread to fill the whole chromosphere. 
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4. Although fast-to-Alfven mode conversion identically vanishes for exactly vertical magnetic field {0 = 0) it can 
be near-total for small 6 and small horizontal wavenumber k, typically for (j) w 90° . 

5. Nevertheless, warm-plasma effects conspire to reflect these Alfven waves back downwards at small 6. 

6. At larger field inclinations {9 > 15°) conversion can still be significant, though at larger k. These Alfven waves 

arc not reflected by warm-plasma effects. 

The role of this study has been to illuminate the basic processes and dependences on wavenumber, magnetic field 
inclination, and wave attack direction of fast-to-Alfven mode conversion. Subsequent work will build these insights 
into more elaborate numerical simulations of sunspots. Currently, all seismic and convective numerical sunspot models 
unrealistically limit the Alfven speed in the overlying atmosphere for numerical reasons {e.g., Hanasogc 2008; Cameron, 
Gizon, & Duvall 2008; Rempel, Schiissler, & Knolkcr 2009). The results presented here have implications for such 
models, as well as for our understanding of atmospheric waves. In particular, we postulate that the returning fast wave 
is important in understanding seismic travel time data, and therefore needs to me modelled correctly. The conversion 
process may also contribute to the recently-observed coronal transverse oscillations. In this scenario, the full width of 
active region chromospheres may be thought of as source region for these coronal waves, rather than (or in addition 
to) direct photospheric motions. 

Finally, it is of interest to rcficct on the recent modelling of coupling between kink and Alfven waves in simple loop 
models with transverse Alfven speed gradients by Pascoe et al. (2010, 2011). A simple straight 'loop' consisting of 
a uniform core and linear decreasing-density shoulders is shaken at its footpoint setting up an 'Alfven' wave in the 
core that rapidly decays by resonant absorption to leave oscillations restricted to the shoulders (see Pascoe et al. 2010, 
Figure 5). However, such back-and-forth footpoint shaking preferentially initiates m = 1 disturbances in the cylinder. 
Strictly, only m = (torsional) motions are genuine incompressive Alfven waves in the cylindrical context, and these do 
not leak into the resonance as they cannot transport energy across field lines. Rather, the rn = 1 'Alfvcnic' oscillations 
in the cores are fast waves, and their absorption at the resonance is in essence the process discussed in Cally & Andries 
(2010) for the case of magnetic field transverse to the direction of inhomogeneity. The fact that Pascoe et al. see 
essentially total absorption rather than the 50% maximum found by Cally & Andries results from multiple absorptions 
as the fast wave bounces back-and-forth across the tube. This is recognized in the schematic Figure 2 of Pascoe et 
al. (2011). In a zero /3 plasma the core fast waves are the so-called (but unfortunately named) compressional Alfven 
waves (Priest 1982, Section 4.3.2). 

APPENDIX 

DETAILS OF NUMERICAL SOLUTION METHOD 

Analytic expressions for the fast and Alfven solutions are required as a; ^ ±oo in order to provide boundary 
conditions for numerical solution of the differential equations and to disentangle the two fast and two Alfven modes. 

WKB at Large s 

A WKB solution is appropriate on ,s = e~^^'^ ^ 1 + since there k^h ^ 1. Substitute the eikonal ansatz 
U = e**^ Yl^=o ^^^^ ~ "^U, assuming Vq ^ Vi ^ V2 ^ . . . and each varying slowly compared to ip. To 
lowest order AVq = is(fi'Yo. Hence isip' and Vg are respectively the eigenvalues and eigenvectors of A. Equivalently, 
the a;-wavenumbers Kx = —s(p' are the eigenvalues of iA, Kx = hdip/dx = —sip' = \/s — k'^, —\/s — k^, s^/"^ secO — 
Kiaa.6 cos(j), and — s^/^ sec 9 — Ktan0cos</), representing respectively the right-moving fast wave, the left-moving fast 
wave, the rightward Alfven wave, and the leftward Alfven wave. Of course, Vq is determined by this process only up 
to a multiplicative scalar function of s. To find this scalar we must progress to the next order. 

Adapting and completing the method of Weinberg (1962) (see his Equation (131)), we have 

(A-is<^'I)V„ = sV;_i. (Al) 

The amplitude dependence of Vq on s may then be determined by setting Vg = /(s)vo, where vq is an arbitrarily 
scaled eigenvector of A, and then premultiplying Equation (Al), for n = 1, by w^, where w is the corresponding left 
eigenvector, thereby extinguishing the left hand side. This leaves (In/)' = --w-^Vq/w-^Vq which may be integrated 
exactly to find /(s). 

The V„ (n ^ 1) may then be calculated recursively by solving Equation (Al) with constraint w^V^ = 0. Let v„ 
be a particular solution of Equation (Al), with the general solution being V„ = v„ -|- a„(s)Vo. Since V„ <i; Vq is 
required as s — > 00, the constraint fixes q:„ = J"" W"^v', /w"^Vo ds', with the integrand assumed a function of the 
dummy variable s' . This integral must be evaluated numerically in practice. Excellent results are attained with the 
expansion truncated at n = 1. It is found that Vi = 0(s~^/^)Vo as s — >■ 00, so the Vi <C! Vq assumption is indeed 
met on s 1. 

Frobenius about s = 

Ks X ^ +00 we utilize the Frobenius expansion U = X^^q about s = 0. This has an infinite radius of 

convergence since there are no singularities other than s = itself. Substituting this expansion into Equation (5) yields 
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the eigenvalue indicial equation AqUq = /iUo, thereby specifying both the hitherto unknown index /i and the zeroth 
coefRcient uq as the eigenvalues and eigenvectors of Ao respectively. Later coefficients are found by the recursion 
u„ = — (Ao — (n + /Li)I)~^AiU„_i. The eigenvalues are /U = — k, k, and i k tan cos </>, representing respectively (as 
s — )• 0+) the exponentially growing fast wave (Ui), the decaying fast wave (U2), and an Alfven wave (U3). Specifically, 

for the physical fast eigenvalue /x = /i2 = the eigenvector is = {— cos cos (j) + i sin 6, sin (j),—K cos 9 cos (j) + 
iKsinO, Ksin <p}. 

The Alfvenic eigenvalue /is has algebraic multiplicity 2 but geometric multiplicity 1. The only independent eigen- 

vector is Uq = {— i cos cot tan 0, — i cot 6*, kcos sin kcos (^}. The series expansion has failed to find the full 
complement of four independent solutions. The fourth solution must therefore take the form 

00 

U4 = U3 In s + ^ v„ 5"+"^ , (A2) 

n=0 

where (Aq — /X3l)vo = Uq and 

(Ao - (n + /i3)I)v„ = (u„ - Aiv„_i) , (n>0). (A3) 

Here the u„ vectors are the coefficients appearing in U3. The equation for vo is of course singular and admits an 
arbitrary additive multiple of Uq in its solution. The later v„ are fully determined by Equation (A3) in general. 

It remains to determine the combination of U3 and U4 (or equivalcntly the additive multiple of U3 required in U4) 
that corresponds to the pure outgoing Alfven wave. Note that Equation (1) shows that — VpX is the source term for 
otherwise free Alfven waves propagating along field lines. As shown by Cally & Andries (2010), this driving of Alfven 
waves occurs in a compact mode conversion region near the reflection point, and is characterized by a stationary phase 
integral. Well to the right of this (s k) x is both very small and of the wrong phase to interact. To a good ap- 
proximation therefore, the solution in s <C k is just that of + uj^/a^^ € = 0, viz., ^ oc s^'^*'^'^^'^°^'^Hq^''^\2^sscc6), 

where H^^^ is the leftward propagating Hankel function (increasing s) and H^'^^ moves to the right. We retain the 
latter only. Now 7?^^^ (2^5 sec 0) = 1 -i[2(lnsec 6' + C) + ln .s]/7r + e'(.s), where C = 0.577216 ... is Euler's constant. This 
is sufficient to determine the radiating Alfven solution, U5 = KV^ + U4, where K = in + Inscc^ 6 + 2C — d-vg/d-Uo- 
Here, the tilde denotes the restriction of a four- vector to its first two components only, and d = (cos 9 sin </>, cos </>) (in 
(_L,j/) space) is the polarization direction of U4 as s ^ 0. The idea is to project U3 and U4 onto their asymptotic 
polarization direction, and force the coefficients of 1 and Ins to appear in the same proportion as in the Hankel 
function. It may be confirmed that V- U5 = to leading order in the Probenius expansions (n = 0), but not beyond 
(X5 = 0{s^^'^^ Ins)), indicating that U5 is a 'pure' Alfven wave only in the limit s ^ 0. 

Shooting Method 

A bi-directional shooting method is used to solve Equation (5) subject to the above boundary conditions. Because 
of the very different eigenvalues of the four solutions, an automatically switching stiff/nonstiff method is used. Despite 
this, resolution is lost and the system becomes ill-conditioned as k increases, typically beyond about 8. This limit can 
be pushed a little higher using high precision arithmetic (20 - 40 decimal digits), but this soon becomes prohibitively 
expensive. Thankfully, waves of interest in the solar atmosphere do not usually have wavelengths small compared to 
the density scale height, so «; 1 is of little practical concern. 

The authors thank Scott Mcintosh for a careful reading of the manuscript and valuable suggestions. 
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SUPPLEMENTARY MATERIAL 
Tables 2 to 8 are to be included in 'Supplementary Material', not in the paper itself. 
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Table 2 

Forward Alfven conversion coefficient J?+ for a range of k and (j> at magnetic field inclination 6 = 20°. 
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,106 





.265 


0.445 


0, 


,596 


0, 


,681 


0, 


,685 


0, 


,620 


0, 


,510 


0, 


,386 


0, 


,271 


0, 


,177 


0, 


,108 


0, 


,062 


0, 


,032 


0.8 


0.013 


0, 


,111 





.277 


0.459 


0, 


,605 


0, 


,676 


0, 


,663 


0, 


,581 


0, 


,462 


0, 


,336 


0, 


,226 





,141 





,083 


0, 


,045 


0, 


,023 


0.9 


0.014 


0, 


,116 





.288 


0.471 


0, 


,611 


0, 


,670 


0, 


,640 


0, 


,544 


0, 


,417 


0, 


,291 


0, 


,187 





,112 





,063 


0, 


,033 


0, 


,016 


1.0 


0.014 


0, 


,121 





.298 


0.483 


0, 


,616 


0, 


,662 


0, 


,616 


0, 


,508 


0, 


,375 


0, 


,252 


0, 


,155 


0, 


,089 


0, 


,048 


0, 


,024 


0, 


,011 


1.5 


0.017 


0, 


,144 


0, 


.341 


0.525 


0, 


.623 


0, 


,605 


0, 


,497 


0, 


,352 


0, 


,218 


0, 


,120 


0, 


,060 


0, 


.027 


0, 


.012 


0, 


,005 


0, 


,002 


2.0 


0.020 


0, 


,163 





.376 


0.551 


0, 


.608 


0, 


,538 


0, 


,391 


0, 


,238 


0, 


,123 


0, 


,056 


0, 


,022 


0, 


.008 


0, 


.003 


0, 


,001 


0, 


,000 


2.5 


0.022 


0. 


,181 





.404 


0.566 


0, 


.583 


0, 


.470 


0, 


.303 


0, 


,159 


0, 


,069 


0, 


,025 


0, 


,008 


0, 


.002 


0, 


.001 


0, 


.000 


0, 


.000 


3.0 


0.025 


0. 


,196 





.428 


0.572 


0, 


.552 


0, 


.406 


0, 


.232 


0, 


,105 


0, 


,038 


0, 


,012 


0, 


,003 


0, 


.001 


0, 


.000 


0, 


.000 


0, 


.000 


3.5 


0.027 


0. 


,211 





.448 


0.573 


0, 


.518 


0, 


.348 


0, 


.177 


0, 


,069 


0, 


,021 


0, 


,005 


0, 


,001 


0, 


.000 


0, 


.000 


0, 


.000 


0, 


.000 


4.0 


0.029 


0. 


,224 





.464 


0.570 


0, 


.483 


0, 


.297 


0, 


.134 


0, 


,045 


0, 


,012 


0, 


,002 


0, 


,000 


0, 


.000 


0, 


.000 


0, 


.000 


0, 


.000 


4.5 


0.031 


0. 


,237 





.478 


0.563 





.448 


0, 


.252 


0, 


.101 


0, 


,030 


0, 


,006 


0, 


,001 


0, 


,000 


0, 


.000 





.000 


0, 


.000 


0, 


.000 


5.0 


0.033 


0. 


,249 





.490 


0.554 





.414 


0, 


.213 


0, 


.076 


0, 


,019 


0, 


,003 


0, 


,000 


0, 


,000 


0, 


.000 





.000 


0, 


.000 


0, 


.000 



Table 3 

Forward Alfven conversion coefficient J?+ for a range of k and 4> at magnetic field inclination 9 = 30° . 

















Polarization angle 4> 






K 


5° 


15° 


25° 


35° 


45° 


55° 


65° 75° 85° 


95° 


105° 115° 125° 135° 145° 



0.000 


0.000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 





.000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


0.001 


0.005 


0, 


,010 


0, 


,004 


0, 


,002 


0, 


,001 


0, 


,001 


0, 


,001 


0, 


,001 


0, 


,000 


0, 


,000 





.001 


0, 


,001 


0, 


,001 


0, 


,001 


0, 


,002 


0.01 


0.006 


0, 


,049 


0, 


,105 


0, 


,118 


0, 


,099 


0, 


,076 


0, 


,060 


0, 


,050 


0, 


,046 


0, 


,045 





.049 


0, 


,058 


0, 


,073 


0, 


,094 


0, 


,111 


0.1 


0.007 


0, 


,060 


0, 


,158 


0, 


,283 


0, 


,411 


0, 


,520 


0, 


,591 


0, 


,619 


0, 


,611 


0, 


,573 





,513 





,435 


0, 


,343 


0, 


,246 


0, 


,156 


0.2 


0.008 


0, 


,070 


0, 


,181 


0, 


,319 


0, 


,456 


0, 


,567 


0, 


,631 


0, 


,643 


0, 


,607 


0, 


,535 





.442 





,342 





,247 


0, 


,164 


0, 


,097 


0.3 


0.009 


0, 


,079 


0, 


,201 


0, 


,349 


0, 


,487 


0, 


,587 


0, 


,630 


0, 


,613 


0, 


,547 


0, 


,453 





,349 





,251 





,168 


0, 


,105 


0, 


,059 


0.4 


0.010 


0, 


,087 


0, 


,220 





,375 





,511 


0, 


,597 


0, 


,615 


0, 


,571 


0, 


,482 


0, 


,374 





.269 





,181 


0, 


,113 


0, 


,066 


0, 


,035 


0.5 


0.011 


0, 


,095 


0, 


,238 


0, 


,398 


0, 


,530 


0, 


,599 


0, 


,594 


0, 


,525 


0, 


,420 


0, 


,306 


0, 


.205 


0, 


,128 


0, 


,075 


0, 


,041 


0, 


,020 


0.6 


0.012 


0, 


,103 


0, 


,254 


0, 


.418 


0, 


.544 


0, 


.596 


0, 


.568 


0, 


,479 


0, 


,362 


0, 


,248 





.155 


0, 


.090 


0, 


.049 


0, 


.025 


0, 


.012 


0.7 


0.013 


0, 


,110 


0, 


,269 


0, 


.436 


0, 


.554 


0, 


.589 


0, 


.540 


0, 


,434 


0, 


,310 


0, 


,199 





.117 


0, 


.063 


0, 


.032 


0, 


.015 


0, 


.007 


0.8 


0.014 


0, 


,117 


0, 


,283 


0, 


.451 


0, 


.561 


0, 


.579 


0, 


.511 


0, 


,392 


0, 


,265 


0, 


,160 





.087 


0, 


.044 


0, 


.021 


0, 


.009 


0, 


.004 


0.9 


0.015 


0, 


,123 


0, 


,296 





.465 





.565 


0, 


.567 


0, 


.481 


0, 


,352 


0, 


,225 


0, 


,127 





.065 





.030 





.013 


0, 


,006 


0, 


,002 


1.0 


0.016 


0, 


,129 


0, 


,309 


0, 


,477 


0, 


,568 


0, 


,552 


0, 


,452 


0, 


,315 


0, 


,191 


0, 


,101 





.048 


0, 


,021 


0, 


,009 


0, 


,003 


0, 


,001 


1.5 


0.019 


0, 


,158 


0, 


,360 


0, 


,518 


0, 


,555 


0, 


,468 


0, 


,318 


0, 


,176 


0, 


,081 


0, 


,031 





.011 


0, 


,003 


0, 


,001 


0, 


,000 


0, 


,000 


2.0 


0.023 


0, 


,182 


0, 


,399 


0, 


,537 


0, 


,521 


0, 


,382 


0, 


,216 


0, 


,095 


0, 


,033 


0, 


,009 





.002 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


2.5 


0.026 


0, 


,203 


0, 


,429 


0, 


,541 


0, 


,477 


0, 


,306 


0, 


,144 


0, 


,051 


0, 


,013 


0, 


,003 





.000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


3.0 


0.029 


0, 


,223 


0, 


,452 


0, 


,537 


0, 


,430 


0, 


,241 


0, 


,095 


0, 


,027 


0, 


,005 


0, 


,001 





.000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


3.5 


0.031 


0, 


,240 


0, 


,471 


0, 


,526 


0, 


,384 


0, 


,189 


0, 


,062 


0, 


,014 


0, 


,002 


0, 


,000 





.000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


4.0 


0.034 


0, 


,256 


0, 


,486 


0, 


,511 


0, 


,341 


0, 


,147 


0, 


,041 


0, 


,007 


0, 


,001 


0, 


,000 





.000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


4.5 


0.037 


0, 


,271 


0, 


,497 





,494 





,301 


0, 


,114 


0, 


,026 


0, 


,004 


0, 


,000 


0, 


,000 





.000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


5.0 


0.039 


0, 


,284 


0, 


,505 


0, 


,474 


0, 


,264 


0, 


,088 


0, 


,017 


0, 


,002 


0, 


,000 


0, 


,000 


0, 


.000 


0, 


.000 


0, 


.000 


0, 


,000 


0, 


,000 
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Table 4 

Forward Alfven conversion coefficient J?+ for a range of k and (j> at magnetic field inclination 6 = 40°. 

















Polarization angle (j) 












K 


5° 


15° 


25° 


35° 


45° 


55° 


65° 75° 85° 


95° 


105° 


115° 


125° 135° 


145° 



0. 


.000 


0.000 


0. 


.000 


0.000 


0.000 


0, 


.000 


0.000 


0. 


.000 


0. 


.000 


0. 


.000 


0. 


.000 





.000 


0, 


.000 


0.000 


0.000 


0. 


.000 


0, 


.001 


0.004 


0. 


.005 


0.002 


0.001 


0, 


.000 


0.000 


0. 


.000 


0. 


.000 


0. 


.000 


0. 


.000 





.000 


0, 


.000 


0.000 


0.000 


0. 


.001 


0, 


.01 


0.005 


0. 


.038 


0.070 


0.064 


0, 


.044 


0.029 


0, 


.021 


0. 


.016 


0. 


.014 


0. 


.013 





.015 


0, 


.020 


0.028 


0.041 


0. 


.058 


0, 


.1 


0.006 


0. 


.050 


0.132 


0.233 


0, 


.330 


0.400 


0, 


.427 


0. 


.417 


0. 


.388 


0. 


.354 





.318 


0, 


.274 


0.218 


0.157 


0. 


.098 


0, 


.2 


0.007 


0. 


.061 


0.158 


0.275 





.386 


0.465 


0, 


.496 


0. 


.479 


0. 


.426 


0. 


.355 





.278 


0, 


.204 


0.139 


0.087 


0. 


.049 


0, 


.3 


0.008 


0. 


.072 


0.182 


0.310 





.422 


0.491 


0, 


.501 


0, 


.457 


0. 


.378 


0. 


.287 





.202 





.132 


0.080 


0.045 


0, 


.023 


0, 


.4 


0.010 


0. 


.082 


0.204 


0.340 


0, 


.450 


0.502 


0, 


.487 


0, 


.417 


0, 


.319 


0, 


.221 





.140 


0, 


.082 


0.045 


0.023 


0, 


.011 


0, 


.5 


0.011 


0. 


.091 


0.224 


0.367 


0, 


.470 


0.504 


0, 


.464 


0, 


.373 


0, 


.264 


0, 


.167 





.095 


0, 


.050 


0.025 


0.011 


0, 


.005 


0, 


.6 


0.012 


0. 


.100 


0.243 


0.389 


0, 


.485 


0.500 


0, 


.437 


0, 


.329 


0, 


.215 


0, 


.124 





.064 


0, 


.030 


0.013 


0.006 


0, 


.002 


0, 


.7 


0.013 


0. 


.108 


0.261 


0.409 


0, 


.495 


0.490 


0, 


.407 


0, 


.287 


0, 


.174 


0, 


.091 





.042 


0, 


.018 


0.007 


0.003 


0, 


.001 


0, 


.8 


0.014 


0. 


.116 


0.277 


0.426 


0, 


.501 


0.478 


0, 


.377 


0, 


.249 


0, 


.139 


0, 


.067 





.028 





.011 


0.004 


0.001 


0, 


.000 


0, 


.9 


0.015 


0. 


.124 


0.292 


0.440 


0, 


.504 


0.462 


0, 


.347 


0. 


.215 


0, 


.111 


0, 


.048 





.018 





.006 


0.002 


0.001 


0, 


.000 


1, 


.0 


0.016 


0. 


.131 


0.306 


0.453 


0, 


.505 


0.445 


0, 


.318 


0, 


.184 


0, 


.088 


0, 


.035 





.012 


0, 


.004 


0.001 


0.000 


0, 


.000 


1, 


.5 


0.020 


0. 


.164 


0.362 


0.492 


0, 


.480 


0.352 


0, 


.195 


0, 


.082 


0, 


.027 


0. 


.007 


0, 


.001 


0, 


.000 


0.000 


0.000 


0, 


.000 


2, 


.0 


0.024 


0. 


.192 


0.404 


0.504 


0, 


.433 


0.265 


0, 


.115 


0. 


.035 


0. 


.008 


0. 


.001 





.000 


0, 


.000 


0.000 


0.000 


0, 


.000 


2. 


.5 


0.028 


0. 


.217 


0.434 


0.500 


0, 


.380 


0.194 


0. 


.066 


0. 


.015 


0. 


.002 


0. 


.000 





.000 


0, 


.000 


0.000 


0.000 


0. 


.000 


3. 


.0 


0.032 


0. 


.238 


0.457 


0.487 


0, 


.328 


0.140 


0. 


.038 


0. 


.006 


0. 


.001 


0. 


.000 





.000 


0, 


.000 


0.000 


0.000 


0. 


.000 


3. 


.5 


0.035 


0. 


.258 


0.473 


0.468 


0, 


.279 


0.101 


0. 


.021 


0. 


.003 


0. 


.000 


0. 


.000 





.000 


0, 


.000 


0.000 


0.000 


0. 


.000 


4. 


.0 


0.038 


0. 


.275 


0.485 


0.446 


0, 


.236 


0.071 


0. 


.012 


0. 


.001 


0. 


.000 


0. 


.000 





.000 


0, 


.000 


0.000 


0.000 


0. 


.000 


4, 


.5 


0.041 


0. 


.291 


0.493 


0.421 





.199 


0.051 


0, 


.007 


0. 


.000 


0. 


.000 


0. 


.000 





.000 


0, 


.000 


0.000 


0.000 


0. 


.000 


5, 


.0 


0.044 


0. 


.306 


0.498 


0.396 





.166 


0.036 


0, 


.004 


0. 


.000 


0. 


.000 


0. 


.000 





.000 


0, 


.000 


0.000 


0.000 


0, 


.000 



Table 5 

Forward Alfven conversion coefficient J?+ for a range of k and (f) at magnetic field inclination 6 = 50°. 

















Polarization angle 4> 














K 


5° 


15° 


25° 


35° 


45° 


55° 


65° 75° 85° 


95° 


105° 


115° 


125° 


135° 


145° 



0, 


.000 


0.000 


0, 


.000 


0, 


.000 


0, 


.000 


0, 


.000 


0.000 


0, 


.000 


0, 


.000 


0, 


.000 


0.000 





.000 


0, 


.000 


0, 


.000 


0.000 


0, 


.000 


0, 


.001 


0.003 


0. 


.003 


0. 


.001 


0, 


.001 


0, 


.000 


0.000 


0, 


.000 


0, 


.000 


0, 


.000 


0.000 


0, 


.000 


0, 


.000 


0, 


.000 


0.000 


0, 


.001 


0. 


.01 


0.003 


0. 


.027 


0. 


.049 


0, 


.041 


0, 


.026 


0.015 


0, 


.009 


0, 


.006 


0, 


.005 


0.005 





.006 


0, 


.009 


0, 


.014 


0.023 


0, 


.036 


0. 


.1 


0.005 


0. 


.039 


0. 


.102 





.179 


0, 


.249 


0.289 


0, 


.285 


0. 


.250 


0. 


.211 


0.185 





.169 





.151 


0, 


.122 


0.086 


0, 


.052 


0. 


.2 


0.006 


0. 


.051 


0. 


.130 


0, 


.224 


0, 


.309 


0.360 


0. 


.364 


0. 


.326 


0. 


.266 


0.205 





.150 


0, 


.103 


0, 


.065 


0.037 


0. 


.019 


0. 


.3 


0.007 


0. 


.063 


0. 


.157 


0, 


.264 


0, 


.350 


0.391 


0. 


.376 


0. 


.316 


0. 


.236 


0.159 





.099 


0, 


.056 


0, 


.030 


0.015 


0. 


.007 


0. 


.4 


0.009 


0. 


.074 


0. 


.182 





.298 





.382 


0.406 


0, 


.367 


0. 


.285 


0. 


.192 


0.114 





.060 





.029 





.013 


0.006 


0, 


.002 


0, 


.5 


0.010 


0, 


.084 


0, 


.205 


0, 


.327 


0, 


.405 


0.410 


0, 


.348 


0, 


.248 


0, 


.150 


0.078 





.036 


0, 


.015 


0, 


.006 


0.002 


0, 


.001 


0, 


.6 


0.011 


0, 


.094 


0, 


.226 


0, 


.353 


0, 


.421 


0.407 


0, 


.323 


0, 


.212 


0, 


.115 


0.053 





.021 


0, 


.007 


0, 


.002 


0.001 


0, 


.000 


0, 


.7 


0.012 


0, 


.104 


0, 


.246 


0, 


.374 


0, 


.432 


0.398 


0, 


.296 


0, 


.179 


0, 


.087 


0.035 





.012 


0, 


.004 


0, 


.001 


0.000 


0, 


.000 


0, 


.8 


0.014 


0, 


.113 


0, 


.263 


0, 


.393 


0, 


.439 


0.385 


0, 


.269 


0. 


.149 


0, 


.065 


0.023 





.007 


0, 


.002 


0, 


.000 


0.000 


0, 


.000 


0, 


.9 


0.015 


0, 


.121 


0, 


.280 


0, 


.408 


0, 


.441 


0.370 


0, 


.242 


0, 


.123 


0, 


.049 


0.015 





.004 


0, 


.001 


0, 


.000 


0.000 


0, 


.000 


1, 


.0 


0.016 


0, 


.130 


0, 


.295 


0, 


.421 


0, 


.441 


0.353 


0, 


.217 


0, 


.101 


0, 


.036 


0.010 





.002 


0, 


.000 


0, 


.000 


0.000 


0, 


.000 


1, 


.5 


0.021 


0, 


.166 


0, 


.356 





.459 





.410 


0.262 


0, 


.117 


0, 


.036 


0, 


.007 


0.001 





.000 


0, 


.000 





.000 


0.000 


0, 


.000 


2, 


.0 


0.025 


0, 


,197 


0, 


,399 





.466 





.358 


0.183 


0, 


.060 


0, 


.012 


0, 


.001 


0.000 





.000 





.000 





.000 


0.000 


0, 


.000 


2, 


.5 


0.030 


0, 


.223 


0, 


.429 


0, 


.457 


0, 


.302 


0.124 


0, 


.030 


0, 


.004 


0, 


.000 


0.000 


0, 


.000 


0, 


.000 


0, 


.000 


0.000 


0, 


.000 


3. 


.0 


0.033 


0. 


.247 


0. 


.451 


0, 


.438 


0, 


.250 


0.082 


0. 


.014 


0. 


.001 


0. 


.000 


0.000 





.000 


0, 


.000 


0, 


.000 


0.000 


0. 


.000 


3. 


.5 


0.037 


0. 


.268 


0. 


.465 


0, 


.414 


0, 


.205 


0.054 


0. 


.007 


0. 


.000 


0. 


.000 


0.000 





.000 


0, 


.000 


0, 


.000 


0.000 


0. 


.000 


4. 


.0 


0.040 


0. 


.287 


0. 


.475 


0, 


.386 


0, 


.166 


0.035 


0. 


.003 


0. 


.000 


0. 


.000 


0.000 





.000 


0, 


.000 


0, 


.000 


0.000 


0. 


.000 


4. 


.5 


0.044 


0. 


.304 


0. 


.480 


0, 


.358 


0, 


.133 


0.023 


0. 


.002 


0. 


.000 


0. 


.000 


0.000 





.000 


0, 


.000 


0, 


.000 


0.000 


0. 


.000 


5. 


.0 


0.047 


0. 


.319 


0. 


.482 


0, 


.330 





.107 


0.015 


0, 


.001 


0. 


.000 


0. 


.000 


0.000 





.000 





.000 





.000 


0.000 


0. 


.000 



Fast-to- Alfven Mode Conversion 



Table 6 

Forward Alfven conversion coefficient J?+ for a range of k and (j> at magnetic field inclination 6 = 60°. 

















Polarization angle (j) 












K 


5° 


15° 


25° 


35° 


45° 


55° 


65° 75° 85° 


95° 


105° 


115° 


125° 135° 


145° 



0.000 


0.000 


0. 


.000 


0, 


,000 


0, 


.000 


0, 


.000 


0, 


.000 


0, 


.000 


0, 


,000 


0, 


,000 


0.000 





.000 


0, 


.000 


0.000 


0.000 


0, 


.000 


0.001 


0.002 


0. 


.003 


0, 


,001 


0, 


.000 


0, 


.000 


0, 


.000 


0, 


.000 


0, 


,000 


0, 


,000 


0.000 





.000 


0, 


.000 


0.000 


0.000 


0, 


.000 


0.01 


0.002 


0. 


.017 


0, 


,035 


0, 


.033 


0, 


.020 


0, 


.010 


0, 


.005 


0, 


,003 


0, 


,002 


0.002 





.003 


0, 


.005 


0.009 


0.017 


0, 


.027 


0.1 


0.003 


0. 


.028 


0, 


,073 


0, 


.127 


0, 


.175 


0, 


.199 


0, 


.184 


0, 


,140 


0, 


,100 


0.083 





.079 


0, 


.073 


0.057 


0.038 


0, 


.021 


0.2 


0.005 


0. 


.041 


0, 


,103 


0, 


.174 


0, 


.234 


0, 


.264 


0, 


.252 


0, 


,203 


0, 


,144 


0.098 





.066 


0, 


.040 


0.022 


0.011 


0, 


.005 


0.3 


0.006 


0. 


.053 


0, 


,132 





.217 





.281 


0, 


.300 


0, 


.268 


0, 


,201 


0, 


,128 


0.072 





.037 





.017 


0.007 


0.003 


0, 


,001 


0.4 


0.008 


0, 


,065 


0, 


,159 


0, 


,256 


0, 


,317 


0, 


,319 


0, 


,265 


0, 


,180 


0, 


,100 


0.047 





.019 


0, 


,007 


0.002 


0.001 


0, 


,000 


0.5 


0.009 


0, 


,077 


0, 


,185 


0, 


,288 


0, 


,344 


0, 


,328 


0, 


,252 


0, 


,154 


0, 


,074 


0.029 





.009 


0, 


,003 


0.001 


0.000 


0, 


,000 


0.6 


0.011 


0, 


,088 


0, 


,208 


0, 


,316 


0, 


,363 


0, 


,327 


0, 


,233 


0, 


,128 


0, 


,053 


0.017 





.004 


0, 


,001 


0.000 


0.000 


0, 


,000 


0.7 


0.012 


0, 


,098 


0, 


,229 


0, 


,340 


0, 


,376 


0, 


,321 


0, 


,211 


0, 


,104 


0, 


,038 


0.010 





.002 


0, 


,000 


0.000 


0.000 


0, 


,000 


0.8 


0.013 


0, 


,108 


0, 


,248 


0, 


,360 


0, 


,383 


0, 


,310 


0, 


,189 


0, 


,084 


0, 


,026 


0.006 





.001 


0, 


,000 


0.000 


0.000 


0, 


,000 


0.9 


0.014 


0, 


,117 


0, 


,266 


0, 


,376 


0, 


,386 


0, 


,296 


0, 


,167 


0, 


,066 


0, 


,018 


0.003 





.000 


0, 


,000 


0.000 


0.000 


0, 


,000 


1.0 


0.016 


0, 


,126 


0, 


,282 


0, 


,390 


0, 


,386 


0, 


,280 


0, 


,146 


0, 


,052 


0, 


,012 


0.002 


0, 


.000 


0, 


,000 


0.000 


0.000 


0, 


,000 


1.5 


0.021 


0, 


,165 


0, 


,347 


0, 


.427 


0, 


.354 


0, 


,197 


0, 


,070 


0, 


,015 


0, 


,002 


0.000 


0, 


.000 


0, 


.000 


0.000 


0.000 


0, 


,000 


2.0 


0.026 


0, 


,198 


0, 


,390 


0, 


.432 


0, 


.300 


0, 


,128 


0, 


,031 


0, 


,004 


0, 


,000 


0.000 





.000 


0, 


.000 


0.000 


0.000 


0, 


,000 


2.5 


0.030 


0. 


,227 


0, 


,421 


0, 


.419 


0, 


.245 


0, 


.081 


0, 


.013 


0, 


,001 


0, 


,000 


0.000 





.000 


0, 


.000 


0.000 


0.000 


0, 


.000 


3.0 


0.034 


0. 


,251 


0, 


,441 


0, 


.396 


0, 


.196 


0, 


.050 


0, 


.006 


0, 


,000 


0, 


,000 


0.000 





.000 


0, 


.000 


0.000 


0.000 


0, 


.000 


3.5 


0.038 


0. 


,273 


0, 


,454 


0, 


.369 


0, 


.154 


0, 


.030 


0, 


.002 


0, 


,000 


0, 


,000 


0.000 





.000 


0, 


.000 


0.000 


0.000 


0, 


.000 


4.0 


0.042 


0. 


,293 


0, 


,461 


0, 


.339 


0, 


.121 


0, 


.018 


0, 


.001 


0, 


,000 


0, 


,000 


0.000 





.000 


0, 


.000 


0.000 


0.000 


0, 


.000 


4.5 


0.046 


0. 


,311 


0, 


,465 


0, 


.310 





.093 


0, 


.011 


0, 


.000 


0, 


,000 


0, 


,000 


0.000 





.000 


0, 


.000 


0.000 


0.000 


0, 


.000 



Table 7 

Forward Alfven conversion coefficient J?+ for a range of k and (f) at magnetic field inclination 6 = 70°. 

















Polarization an; 
















K 


5° 


15° 


25° 


35° 


45° 


55° 


65° 75° 


85° 


95° 


105° 


115° 


125° 


135° 


145° 



0.000 


0.000 


0, 


,000 


0, 


,000 


0, 


.000 


0, 


.000 


0, 


,000 


0, 


,000 


0, 


,000 


0.000 


0.000 


0, 


.000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


0.001 


0.001 


0, 


,003 


0, 


,002 


0, 


.001 


0, 


.000 


0, 


,000 


0, 


,000 


0, 


,000 


0.000 


0.000 


0, 


.000 


0, 


.000 


0, 


.000 


0, 


,000 


0, 


,001 


0.01 


0.001 


0, 


,009 


0, 


,021 





.027 


0, 


.020 


0, 


.010 


0, 


.004 


0, 


,001 


0.001 


0.001 





.001 





.004 


0, 


.008 


0, 


.016 


0, 


.020 


0.1 


0.002 


0, 


,018 


0, 


,047 


0, 


.081 


0, 


.111 


0, 


.127 


0, 


.116 


0, 


,077 


0.041 


0.030 





.031 


0, 


.027 


0, 


.018 


0, 


.010 


0, 


.005 


0.2 


0.004 


0, 


,032 


0, 


,079 


0, 


.131 


0, 


.171 


0, 


.184 


0, 


.164 


0, 


,116 


0.064 


0.035 





.019 


0, 


.009 


0, 


.004 


0, 


.001 


0, 


.000 


0.3 


0.005 


0, 


,045 


0, 


111 





.179 





.223 


0, 


.225 


0, 


.184 


0, 


,117 


0.057 


0.023 





.008 





.002 





.001 


0, 


,000 


0, 


,000 


0.4 


0.007 


0, 


,058 


0, 


,141 


0, 


,221 


0, 


,264 


0, 


,251 


0, 


,188 


0, 


,106 


0.043 


0.013 





.003 


0, 


,001 


0, 


,000 


0, 


,000 


0, 


,000 


0.5 


0.008 


0, 


,071 


0, 


,168 


0, 


,257 


0, 


,295 


0, 


,264 


0, 


,181 


0, 


,089 


0.030 


0.007 





.001 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


0.6 


0.010 


0, 


,082 


0, 


,193 


0, 


,287 


0, 


,317 


0, 


,268 


0, 


,168 


0, 


,073 


0.020 


0.003 





.000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


0.7 


0.011 


0, 


,093 


0, 


,215 


0, 


,312 


0, 


,332 


0, 


,264 


0, 


,152 


0, 


,058 


0.013 


0.002 





.000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


0.8 


0.013 


0, 


,103 


0, 


,235 


0, 


,333 


0, 


,340 


0, 


,255 


0, 


,134 


0, 


,045 


0.008 


0.001 





.000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


0.9 


0.014 


0, 


,113 


0, 


,254 


0, 


,351 


0, 


,344 


0, 


,243 


0, 


,117 


0, 


,034 


0.005 


0.000 





.000 


0, 


,000 


0, 


,000 


0, 


,000 


0, 


,000 


1.0 


0.015 


0, 


,123 


0, 


,271 





,365 





,344 


0, 


,229 


0, 


,101 


0, 


,026 


0.003 


0.000 





.000 


0, 


,000 





,000 


0, 


,000 


0, 


,000 


1.5 


0.021 


0, 


,164 


0, 


,337 





,402 





,312 


0, 


,154 


0, 


,043 


0, 


,006 


0.000 


0.000 





.000 





,000 





,000 


0, 


,000 


0, 


,000 


2.0 


0.026 


0, 


,198 


0, 


,382 


0, 


,405 


0, 


,260 


0, 


,095 


0, 


,017 


0, 


,001 


0.000 


0.000 


0, 


.000 


0, 


.000 


0, 


,000 


0, 


,000 


0, 


,000 


2.5 


0.031 


0, 


,227 


0, 


,412 


0, 


.390 


0, 


.207 


0, 


.056 


0, 


.006 


0, 


,000 


0.000 


0.000 





.000 


0, 


.000 


0, 


.000 


0, 


.000 


0, 


.000 


3.0 


0.035 


0, 


,253 


0, 


,431 


0, 


.365 


0, 


.161 


0, 


.033 


0, 


.002 


0, 


,000 


0.000 


0.000 





.000 


0, 


.000 


0, 


.000 


0, 


.000 


0, 


.000 


3.5 


0.039 


0, 


,276 


0, 


,443 





.336 





.123 


0, 


.019 


0, 


.001 


0, 


,000 


0.000 


0.000 





.000 





.000 





.000 


0, 


.000 


0, 


.000 


4.0 


0.043 


0, 


,296 


0, 


,449 





.306 


0, 


.093 


0, 


.011 


0, 


.000 


0, 


,000 


0.000 


0.000 





.000 





.000 


0, 


.000 


0, 


.000 


0, 


.000 



P.S. Cally and S.C. Hansen 



Table 8 

Forward Alfven conversion coefficient J?+ for a range of k and (j> at magnetic field inclination 6 = 80°. 

















Polarization angle (j) 












K 


5° 


15° 


25° 


35° 


45° 


55° 


65° 75° 85° 


95° 


105° 


115° 


125° 135° 


145° 



0. 


.000 


0.000 


0. 


.000 


0. 


.000 


0.000 


0, 


.000 


0.000 


0. 


.000 


0.000 


0.000 


0.000 





.000 


0.000 


0.000 


0.000 


0.000 


0, 


.001 


0.000 


0. 


.002 


0. 


.003 


0.002 


0, 


.001 


0.000 


0. 


.000 


0.000 


0.000 


0.000 





.000 


0.000 


0.000 


0.001 


0.002 


0, 


.01 


0.000 


0. 


.003 


0. 


.007 


0.012 


0, 


.015 


0.013 


0. 


.006 


0.001 


0.000 


0.000 





.001 


0.005 


0.009 


0.009 


0.007 


0, 


.1 


0.001 


0. 


.012 


0. 


.029 


0.049 


0, 


.063 


0.068 


0. 


.061 


0.041 


0.013 


0.007 





.006 


0.003 


0.001 


0.000 


0.000 


0, 


.2 


0.003 


0. 


.026 


0. 


.064 


0.103 


0, 


.129 


0.129 


0. 


.102 


0.059 


0.020 


0.006 





.001 


0.000 


0.000 


0.000 


0.000 


0, 


.3 


0.005 


0. 


.040 


0. 


.098 


0.155 





.188 


0.178 


0, 


.129 


0.064 


0.018 


0.003 





.000 


0.000 


0.000 


0.000 


0.000 


0, 


.4 


0.006 


0. 


.054 


0, 


.129 


0.200 


0, 


.233 


0.210 


0, 


.140 


0.061 


0.013 


0.001 





.000 


0.000 


0.000 


0.000 


0.000 


0, 


.5 


0.008 


0. 


.067 


0, 


.158 


0.238 


0, 


.267 


0.227 


0, 


.139 


0.052 


0.008 


0.000 





.000 


0.000 


0.000 


0.000 


0.000 


0, 


.6 


0.010 


0. 


.079 


0, 


.183 


0.269 


0, 


.290 


0.233 


0, 


.131 


0.042 


0.005 


0.000 





.000 


0.000 


0.000 


0.000 


0.000 


0, 


.7 


0.011 


0. 


.090 


0, 


.206 


0.295 


0, 


.306 


0.231 


0, 


.118 


0.033 


0.003 


0.000 





.000 


0.000 


0.000 


0.000 


0.000 


0, 


.8 


0.012 


0. 


.101 


0, 


.227 


0.317 


0, 


.315 


0.223 


0, 


.104 


0.025 


0.002 


0.000 





.000 


0.000 


0.000 


0.000 


0.000 


0, 


.9 


0.014 


0. 


.111 


0, 


.246 


0.335 


0, 


.319 


0.212 


0, 


.090 


0.018 


0.001 


0.000 





.000 


0.000 


0.000 


0.000 


0.000 


1, 


.0 


0.015 


0. 


.120 


0. 


.263 


0.349 


0, 


.319 


0.200 


0, 


.076 


0.013 


0.000 


0.000 





.000 


0.000 


0.000 


0.000 


0.000 


1, 


.5 


0.021 


0. 


.162 


0. 


.331 


0.386 


0, 


.288 


0.130 


0, 


.030 


0.002 


0.000 


0.000 


0, 


.000 


0.000 


0.000 


0.000 


0.000 


2, 


.0 


0.026 


0. 


.198 


0. 


.376 


0.388 


0, 


.236 


0.078 


0, 


.011 


0.000 


0.000 


0.000 


0, 


.000 


0.000 


0.000 


0.000 


0.000 


2. 


.5 


0.031 


0. 


.228 


0. 


.406 


0.372 


0, 


.185 


0.044 


0. 


.004 


0.000 


0.000 


0.000 





.000 


0.000 


0.000 


0.000 


0.000 


3. 


.0 


0.035 


0. 


.254 


0. 


.425 


0.346 





.141 


0.024 


0, 


.001 


0.000 


0.000 


0.000 





.000 


0.000 


0.000 


0.000 


0.000 


3, 


.5 


0.040 


0. 


.277 


0. 


.436 


0.317 





.106 


0.013 


0, 


.000 


0.000 


0.000 


0.000 





.000 


0.000 


0.000 


0.000 


0.000 



